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characteristic of classicalization. The classicalization radius sets the scale for the onset of 
significant deformations of the collapsing configuration, which result in the formation of 
shock fronts. A characteristic observable feature of the classicalization process is the cre- 
ation of an outgoing field configuration that extends far beyond the classicalization radius. 
This feature develops before the deformed wavepacket reaches distances of the order of the 
fundamental scale. We find that in some models the scattering problem may not have real 
solutions over the whole space at late times. We determine the origin of this behavior and 
discuss the consistency of the underlying models. 
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1 Introduction 

The physics of classicalization concerns the nature of high-energy scattering in certain 
classes of nonrenormalizable field theories. The proposal of refs. [1-4] is that scatter- 
ing can take place at length scales much larger than the typical scale associated with the 
nonrenormalizable terms in the Lagrangian. The behavior is similar to ultra-Planckian 
scattering in gravitational theories, during which a black hole is expected to start forming 
at distances comparable to the Schwarzschild radius. For sufficiently large center-of-mass 
energy, the Schwarzschild radius can be much larger than the Planck length. If the forma- 
tion and subsequent evaporation of a black hole are viewed as parts of a scattering process, 
it becomes clear than the physical length scale determining the cross section can be the 
Schwarzschild radius and not the Planck scale. It has been advocated that similar behavior 
may be observed in some nonrenormalizable field theories [1-4]. In scattering processes, 
the center-of-mass energy introduces a new dimensionful scale in the problem that can be 
used to define the analogue of the Schwarzschild radius. This new length scale is termed 
classicalization radius. It defines the distance at which scattering starts and, for sufficiently 
high center-of-mass energy, can be much larger than the fundamental length scale L* of 
the theory, determined by the coefficient of the non-renormalizable term in the action. The 
scenario has very interesting collider phenomenology [5]. 

An idealized scattering process was proposed in refs. [2, 3] as a means to check the va- 
lidity of the classicalization arguments. The process involves an initial spherical wavepacket 
of very large radius that propagates towards the center of symmetry. The solution of the 
classical equation of motion can provide a clear indication of the onset of scattering. One 
expects that it should be possible to verify that this starts at length scales much larger 
than L*. An approximate analytical solution of the equation of motion was presented in 
refs. [2, 3]. It demonstrates that there are significant corrections to the initial wavepacket 
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at length scales of order r*, with the simultaneous creation of a long tail in the field con- 
figuration, stretching out to large distances. However, an exact solution of the equation 
of motion is still lacking. In ref. [6], it was suggested that the equation of motion can 
be written as a quasi-linear partial differential equation of second order. For the initial 
spherical wavepacket this equation is hyperbolic, supporting the propagation of wavelike 
configurations. However, when the center of the wavepacket approaches r#, the equation 
may become elliptic in certain regions. This change could result in the suppression of wave 
propagation and the emergence of scattering. 

In this work we present numerical solutions of the field equations of motion for two 
theories, described by nonrenormalizable Lagrangians, which are expected to generate be- 
havior characteristic of classicalization. We discuss in detail how the typical features of 
classicalization emerge during the evolution of the idealized spherical wavepacket described 
above. We also point out several novel features that seem to be inherent in scattering pro- 
cesses in the context of such theories. In the following two sections we describe the theories 
in question. In section 4 we present the method we employ for the integration of the equa- 
tions of motion. We also discuss the crucial issue of existence or absence of real solutions 
in certain ranges of the radial coordinate r. In section 5 we present the results of the 
numerical integration. In the final section we discuss the interpretation of our results. 



2 Quartic action 



The prototype scalar theory in which classicalization is expected to occur is described by 
the Lagrangian density 



C= 1 -(d^) 2 -5 1 L 



(2.1) 



4 v v ^ 

We allow for both signs of the higher derivative term by assuming that 5\ can take the 
values 5\ = dbl. The equation of motion of the field <j> is 



d' 1 



1 - 5 1 Li 



0. 



(2.2) 



An idealized scattering process, which displays the characteristic behavior associated 
with classicalization, involves a collapsing spherical wavepacket initially represented by a 
Gaussian of width a centered around a radius tq. It has the form 



<fo(t,r) 



A 



exp 



(r + t- r ) 



(2.3) 



Using perturbation theory, it was shown in refs. [2, 3] that the field configuration is strongly 
deformed when the peak of the wavepacket reaches the classicalization radius 



A 2 L> 



1/3 



(2.4) 



This can be substantially larger than the fundamental length scale L* when the center-of- 
mass energy s ~ A 2 /a is much larger than l/L*. It must be noted that the analysis of refs. 
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[2, 3] was carried out for the case 5± = — 1. However, it is straightforward to check that 
the same arguments apply to the case 5\ = 1 as well. 

Another point of view was presented in ref. [6]. It was argued that the behavior 
associated with classicalization is related to the dynamical change of type of the classical 
equation of motion. For a spherically symmetric configuration, eq. (2.2) takes the form 

(1 - 3Xcfi 2 + Xcfi 2 ) (f> tt - (1 - Xcfi 2 + 3A0 2 ) rr + AX(j> r 4> t (j> tr = ^ (l - Xcfi 2 + Xcfi 2 ) , (2.5) 

where A = 5\L^ and subscripts denote partial derivatives. The above equation can be 
expressed as a conservation law, in the form 

d t [(fit (1 - H\ + \<fi)] - ^d r [r 2 cfi r (1 - Xcfi 2 + Xcfi 2 )] = 0. (2.6) 



It can also be expressed as 



t 



- (cfi 2 1 



+ cfi 2 ) + - {-34 + 4 + Icgcfi 2 ) 



^8 r [r 2 cfi t cfi r (1 - Ht + A0 2 )] = 0. (2.7) 



This last equation is equivalent to the conservation of the energy- momentum tensor result- 
ing from the Lagrangian of eq. (2.1). The quantity 

P = \ + <fir) + \ (-3^1 + 4 + 2$<ft) (2.8) 

is the local energy density. The initial conditions for the solution of the equation of motion 
are of the Cauchy type: 

0(O,r) = to(O,r), d t (b(0,r) = d t M0,r), (2.9) 

where we assume that tq is much larger than any other physical scale. 

The properties of eq. (2.5) become more transparent if it is written in the form of a 
quasilinear second-order partial differential equation: 

A{4>t, <t>r) 4>tt + B{(fi U (fir) (fitr + C((fit, 4>r) 4>rr = T>((fi t , 4>r,f), (2-10) 

with 

A{(fi t ,(fir) = l-3Xcfi 2 + Xcfi 2 (2.11) 

B(cf> t ,(fi r ) = 4X<fit<fi r (2-12) 

C(<fi t ,<fi r ) = - (l - Xcfi 2 + 3Xcfi 2 ) (2.13) 

V(cfi t ,cfi r ,r) = ^ (1 - Xcfi 2 + Xcfi 2 ) . (2.14) 
The type of this partial differential equation is determined by the discriminant 



A = X -(B 2 - 4 AC) = 3 Q - Xcfi 2 + Xcfi 2 )) (1 - Xcfi 2 + Xcfi 2 ) . (2.15) 

For A > the equation is hyperbolic, for A = parabolic, while for A < elliptic. Hyper- 
bolic equations admit wave-like solutions, while elliptic ones do not support propagating 
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solutions. The observation of ref. [6] is that classicalization may be associated with the 
change of the type of the equation during the evolution of the initial configuration. It can 
be shown [6] that, if A, B, C are evaluated for the configuration (2.3), the discriminant 
(2.15) switches sign in a certain r-interval when the wavepacket moves into the vicinity of 
the classicalization radius (2.4). The sign change occurs for both values of 6±. As a result, 
we expect behavior associated with classicalization in both cases. 

In the following sections we focus on the details of the evolution of an initial field 
configuration given by eq. (2.3). We have managed to obtain only a numerical solution of 
the equation of motion with this initial condition, which we present in section 5. However, 
analytical solutions of the quasilinear second-order partial differential equation (2.5) can 
also be found and they are of interest, even though they do not satisfy the initial condition 
(2.3). 

For (ft = 4>(r), eq. (2.6) gives 

&.(1 + A$) = ^, (2.16) 

with c a constant of integration. If we are interested in a localized configuration of finite 
energy, we must select the root of the above equation for which cf> r approaches c/r 2 for 
large r. The solution with this asymptotic behavior extends down to r = for A > 0. For 
A < the solution displays a square-root singularity at a certain value r s ^ 0. It is possible 
to join smoothly this solution with a second root of eq. (2.16) that again extends from 
r s to infinite r. In this way the field (ft becomes a double-valued function of position. A 
similar construction was employed in ref. [7] in order to describe Dirichlet (d — l)-branes 
embedded in (d + l)-dimensional Minkowski space, with the field (ft corresponding to the 
transverse coordinate. However, it is difficult to find an interepretation of such a solution 
within the framework of the four-dimensional scalar field theory that we are considering: 
The field (ft would be double- valued for r > r s , while it would not be defined for r < r s . 

The partial differential equation (2.5) also has exact dynamical solutions of the form 
(ft = 4>{z), with z = r 2 — (t — to) 2 . They are given by 

(ft t (t, r) = -2(t - t ) h(r 2 - (t - t ) 2 ) (2.17) 
(ft r (t,r) = 2rh(r 2 -(t-to) 2 ), (2.18) 

with h(z) satisfying the differential equation 

z(l + 12A z h(z) 2 )h'(z) + 12A z h{zf + 2h(z) = 0. (2.19) 

This equation has three branches of solutions, similarly to eq. (2.16). They correspond to 
the roots of the cubic equation 

4Xz 2 h 3 (z) + z 2 h(z) + c = 0, (2.20) 

with c an integration constant. We do not present the explicit expressions because of their 
complexity. In the following section, we discuss in detail similar solutions that can be 
obtained in the context of a theory of the Dirac-Born-Infeld type. 
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3 Dirac-Born-Infeld action 



A related theory, which can also display the classicalization phenomenon, is described by 
a Lagrangian density of the Dirac-Born-Infeld (DBI) type: 



C 



1 



1 - 6 2 L* (8^ 



(3.1) 



with ^2 = ±1, as in the previous case. The equation of motion of the field <j> is 



0'' 



0„ 



1 - 5 2 Li (d^ 



(3.2) 



The leading terms in the expansion of the l.h.s. of the above equations reproduce the 
quartic model of the previous section, provided that we take 5 2 = —25±. The Lagrangian 
(3.1) also includes a cosmological constant term, which is irrelevant for our considerations. 
In general, we expect similar behavior for the solutions of eqs. (2.2), (3.2) if <5i, 5 2 are 
taken with opposite signs. However, the higher-order terms resulting from the expansion 
of the l.h.s. of eq. (3.2) are also important for the dynamics. 

When expressed in spherical coordinates, eq. (3.2) can be put in the form 



(1 + A$) <p t 



1 



A0f) 4>r 



2\<p r 4> t 4> tT = ^ (1 - Xtf + Xcfi) , (3.3) 



where A = 5 2 L\ and subscripts denote partial derivatives. It is important to emphasize that, 
in order to obtain the above equation, eq. (3.2) must be multiplied by (l — X(f> 2 + \<p 2 ) 3 ^ 2 . 
As a result, we must assume that 



1 - \$ + \<p 2 . > 0. 



(3.4) 



This assumption is an obvious constraint imposed by the form of the Lagrangian density 
(3.1). Eq. (3.3) can also be expressed as a conservation law, in the form 



dt 



It can also be written in the form 



dt 



~2 dr 



y/1 ~ Ml + H 2 r 



'->r<Pt 



o. 



o. 



(3.5) 



(3.6) 



This last equation is equivalent to the conservation of the energy-momentum tensor result- 
ing from the Lagrangian density of eq. (3.1). The quantity 



P 



A0! 



1 

A 



(3.7) 



is the local energy density, after the subtraction of a cosmological constant included in the 
Lagrangian. 
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Eq. (3.3) can be written in the form (2.10) with 



A{4>tAr) = I + \4? T (3-8) 

B(&,&.) = -2A&^. (3.9) 
C(<k,<f>r) = - (1 - A0f ) (3.10) 
X>(&, 4> r , r) = ^ (1 - A0 2 + X<f> 2 r ) . (3.11) 

The type of this partial differential equation is determined by the discriminant 

A = ^(B 2 -4AC) = 1- X$ + X$. (3.12) 

For A > the equation is hyperbolic, for A = parabolic, while for A < elliptic. The 
DBI theory represents a special case. Requiring the Lagrangian density to be real results 
in the constraint (3.4), which excludes the possibility A < 0. We can again estimate the 
distance r at which eq. (2.10) changes type if A, B, C are evaluated for the configuration 
(2.3). For both values of 82, the discriminant (3.12) may vanish when the wavepacket 
approaches a region in which r is approximately given by eq. (2.4) [6]. 

Exact analytical solutions of the evolution equation can be derived, similarly to the 
previous section. For <\> = cp(r), eq. (3.5) gives <f> r / a/1 + X(j) 2 = ±c/r 2 , where we have 
assumed that the constant of integration c is positive. This relation can be written as 

(3.13) 



Vr 4 - Ac 2 ' 

In ref. [1] this solution was associated with the phenomenon of classicalization. It was 
interpreted as a field configuration induced by a <5-function source resulting from the large 
concentration of energy within a small region of space around r = 0. For A < the solutions 
for both signs extend down to r = 0. For A > they display a square-root singluarity at 
r s = AVW They can be joined smoothly in order to create a continuous double-valued 
function of r that extends from infinite r to r s and back out to infinity. 

Exact dynamical solutions, of the form </> = 4>(z), with z = r 2 — (t — to) 2 , can also be 
found. They are described by eqs. (2.17), (2.18), with h(z) satisfying 

zh'(z) + 6Xzh(zf + 2h(z) = 0. (3.14) 

The solutions of this equation are 

= ±-=L==. (3.15) 

Vcz 4 - AXz 

Requiring that the solution remain real for r — > 00 imposes c > 0. For both signs of A, the 
two solutions display square-root singularities at the value z s = r 2 — (t s — to) 2 that satisfies 
zl = 4A/c. 

An interesting question is whether the analytical solutions we described capture some 
features of the evolution of the incoming spherical wavepacket in the scattering problem. 
Unfortunately, we have not found any concrete evidence for this possibility. In ref. [8], 
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solutions analogous to eqs. (2.17), (2.18), (3.15) were employed in order to describe shock 
waves in the context of lower-dimensional field theories. The shock waves we shall encounter 
in the following sections can be fitted by functions with square-root singularities, but 
with coefficients that do not match the ones deduced from eq. (3.15). It seems that the 
connection of the known static or dynamical analytical solutions with the phenomenon 
of classicalization must be considered rather sketchy. However, it is possible that richer 
analytical solutions may make the connection more concrete. 



4 Numerical method 

The numerical integration of the equations of motion is achieved through the implemen- 
tation of a variant of the leap-frog scheme. We have found that the most efficient form of 
the equations is given by eqs. (2.6), (3.5). The discretized version that we employ is 

Ar 



—r 2 
At 1 



G W + \v j+1 ) -G ut\V/ 



(4.1) 

with U = <k,V = 4> r - The functions F and G are F(U, V) = U{1 - XU 2 + XV 2 ), G{U, V) = 
V(l - XU 2 + XV 2 ) for the first model and F(U, V) = U(l - XU 2 + XV 2 )" 1 / 2 , G(U, V) = 
V(l — XU 2 + XV 2 )~ l l 2 for the second. The index j determines the discretized values of 
the time coordinate t, while the index i the discretized values of the radial coordinate r. 
The above equation must be complemented by an equation that enforces the condition 
d(j)t/dr = d(f> r /dt. Its discretized version is 

Ui +1 -Ui_ 1 = %{yt 1 -Vt')- (4-2) 

In practice, we solve eq. (4.2) for Vf + , which we then substitute in eq. (4.1). The solution 
of the resulting algebraic equation leads to the determination of U- +1 in terms of known 
values at earlier times. The boundary conditions at the ends of the radial interval are not 
of high importance, as our analysis is restricted to time intervals during which the field is 
negligible there. 

The above scheme provides an accurate description of the evolution of the field configu- 
ration, starting from an initial condition given by eq. (2.3) with ro 3> 3> L*. The initial 
stage is simple, as the equation of motion is essentially the wave equation. The evolution 
consists mainly of the displacement of the peak of the wavepacket towards smaller values 
of r, with an increase of its amplitude. The role of the higher-derivative terms becomes 
important during the later stages, when the wavepacket approaches the classicalization 
radius r*. The most problematic aspect for our analysis is that continuous real solutions 
may cease to exist at some stage of the evolution. There are three reasons for this behavior: 

• At some stage the solution develops a shock front. From this point on, the numerical 
integration cannot be continued, as the evolution of the shock depends on additional 
physical assumptions about its nature. The equation of motion becomes singular at 
the location of the shock and the resulting discontinuities in the field configuration, 
or its derivatives, cannot be determined without extra input. We do not attempt 
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Figure 1. The nonlinear wavepacket at various times (solid lines) vs. the linear wavepacket (dotted 
lines), in the context of the DBI theory, described by eq. (3.1), with 82 = 1, = 1. The initial 
wavepacket has A = 20, a = 1. The vertical dashed line denotes the classicalization radius. 



to impose additional conditions, so as to continue the integration, for two reasons: 
Firstly, our intuition on the physical properties of this system is limited; secondly, the 
scattering properties of the solution can be established already before the appearance 
of the shock front. We also mention at this point that the numerical solution is slightly 
better behaved near the shock front, if one employs the quantities F and G, appearing 
in eqs. (4.1), (4.2), as independent variables, instead of U and V. The reformulation 
of the system (4.1), (4.2) is straightforward. Our results have been obtained and 
cross-checked through both formulations. 

• At some time a real solution ceases to exist within a certain range of r. This possibility 
is already apparent in the analytical solutions of the equations of motion that we 
discussed at the ends of sections 2 and 3. We encounter specific examples in the 
numerical analysis in the following section. In the final section we provide an intuitive 
understanding of the origin of this behavior in the context of toy models. 

• The third type of complication occurs at the time when the partial differential equa- 
tion switches type within a range of r. This change is indicated by the discriminant 
A approaching zero from positive values. For A < the equation becomes elliptic 
and its solution requires (Dirichlet or Neumann) boundary conditions on a closed 
contour around the region of interest. Unfortunately, the scattering problem that 
we are considering cannot provide such conditions, as it is set up through Cauchly 
boundary conditions at the initial time. For A < the coordinate t becomes essen- 
tially spatial. Boundary conditions on a closed contour would require the values of 
<j> or its derivatives at times later than the time of interest. In the standard (and 
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2 4 6 8 r 10 12 

Figure 2. The derivatives <j>t (dashed) and <fi r (solid) of the nonlinear field, and the discriminant A 
of eq. (3.12) (solid grey), at two different times, before and after the crossing of the classicalization 
radius. The model is the same as in fig. 1. The vertical dashed line denotes the classicalization 
radius. 

rather limited) studies of partial differential equations of mixed type, the boundary 
conditions are imposed on appropriate boundaries so that a unique solution exists. 
In the case we are considering this is not possible. 

The above complications indicate that the idealized scattering problem that we are 
considering is not always well posed mathematically. The most disturbing issue is the 
inability to find a solution in the range of r where the equation becomes elliptic. A re- 
formulation of the problem is needed in order to overcome this difficulty. Despite these 
misgivings, a solution of the equation of motion over sufficiently long time intervals, so as 
to probe the region of classicalization, is possible in many cases. In the following section 
we present explicit examples. 

5 Numerical results 

We have found that the solutions with the most regular behavior are obtained in the 
context of the DBI model, described by the Lagrangian (3.1), with 82 = 1 and, therefore, 
A > 0. In fig. 1 we depict the solution of the equation of motion at various times. We have 
normalized all dimensonful quantities with respect to the fundamental scale of the theory. 
This is equivalent to setting L* = 1. We have considered an initial spherical wavepacket 
with a width equal to the fundamental scale (a = 1), as the preparation of a narrower 
configuration would presuppose that scales below L* can be probed. The amplitude of the 
configuration must be much larger than 1 in order for the classicalization radius to satisfy 
3> £*. We use A = 20. The solid lines denote the nonlinear field 4>(t, r) resulting from the 
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integration of eq. (2.6). The dotted lines denote the evolution of the initial configuration 
fa{t,r) if the nonlinear terms are neglected (A = 0). In this case the equation of motion is 
the wave equation in spherical coordinates, and eq. (2.3) provides its exact solution. 

As expected, the nonlinear and linear field configurations coincide at early times, when 
the peak of the wavepacket is located at distances much larger than r*. The nonlinear con- 
figuration is subject to visible deformations when the peak approaches r*. The deformations 
become significant when the wavepacket crosses inside the classicalization radius. The final 
configuration in fig. 1 depicts a wavepacket whose front end is located at distances below 
= 1. Its height is reduced significantly in relation to the linear configuration. Its back 
end develops a long tail that extends far beyond the classicalization radius. The last impor- 
tant feature of the configuration is that the part behind the peak is very steep, essentially 
vertical. This indicates the apperance of a shock front, where the field derivatives become 
singular. 

In fig. 2 we display some characteristics of the solution before and after the rvcrossing. 
The dashed and solid lines that almost coincide depict the derivatives fa and fa, respec- 
tively. The grey solid line depicts the discriminant A of eq. (3.3), given by eq. (3.12). The 
form of the field at early times, seen on the right of the vertical dashed line marking the 
classicalization radius, is typical of an incoming spherical wavepacket. The derivatives fa, 
fa are almost equal, while the discriminant deviates only slightly from 1. The late-time 
configuration, located within the classicalization radius, has a significantly different form. 
Within a small r-interval around the shock front, the radial derivative 4> r becomes negative, 
with an absolute value that becomes very large. The time-derivative fa has a similar form. 
At large r, both fa, <j) r approach zero, having opposite signs. This is an important feature 
that we shall discuss in detail below. The discriminant displays three types of behavior: In 
the regions where fa, fa are small, it takes values very close to 1, as the nonlinear terms 
in the equation of motion are subleading. At the location of the shock front, it becomes 
very large. In the remaining regions it approaches zero. This last feature is clearly visible 
in the region 2 ^ r & 3.5. It indicates that the partial differential equation (3.3) becomes 
essentially parabolic within this region. For wavepackets with larger amplitude A, the 
transitions in the values of A between 0, 1 and very large values at the shock front become 
very sharp. 

In fig. 3 we display some additional features of the field configuration before and after 
the rvcrossing. The solid lines depict the field. The late-time configuration is characterized 
by a very long tail that extends well beyond the classicalization radius. The dashed lines 
depict the product 4irr 2 p, with the energy density p given by eq. (3.7). The total energy 
of the configuration, corresponding to the area below the dashed lines, remains constant 
during the evolution. Actually, the conservation of the energy is one of the criteria that we 
employ in order to estimate the accuracy of the numerical integration. For both configu- 
rations depicted in fig. 3, the energy is distributed roughly equally between the front and 
back end of the wavepacket. When the shock front develops, the energy density becomes 
strongly peaked at its location. The long tail of the late-time configuration carries very 
little energy, as both fa, fa are very small there. 

The properties of the tail appearing at late times in figs. 1-3 are the main observable 
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2 4 6 8 10 12 

r 



Figure 3. The nonlinear field cj) (solid) and the product 47rr 2 p, with p the energy density (dashed). 
The model is the same as in fig. 1. The vertical dashed line denotes the classicalization radius. The 
energy density is multiplied by 5 x 1CP 4 . 



consequence of the classicalization scenario. In the region within which the tail develops, 
the equation of motion is essentially the wave equation, as (fit, <j> r are very small and the 
nonlinear corrections negligible. For this reason we expect that the field can be expressed 
as 

<f>{t,r) = -[f(r + t) + g(r-t)}. (5.1) 
r 

The unknown functions / and g can be determined from our numerical solution through 
the relations 

f'(r + t) = I [r(Mt, r) + &(t, r)) + </>(t, r)] (5.2) 

g'{r - t) = - [r{Mt, r) - &(t, r)) + <f>(t, r)] , (5.3) 

where the primes denote derivatives with respect to the arguments. Our results indicate 
that to a very good accuracy f'(r + t) =0, so that we can set consistently f(r + t) = 
in eq. (5.1). Thus the tail of the field represents an outgoing configuration, which can be 
interpreted as a result of the scattering of the initial incoming wavepacket. We have also 
determined the form of the derivative of the function g(r + i), by using data at different 
times. The result is depicted in fig. 4, where we have indicated parts of the curve computed 
at different times. They all form a continuous curve, with overlapping parts that are 
consistent with each to a high accuracy. Within the tail, the field is a monotonically 
increasing function of r+t, starting from negative values and vanishing for large arguments. 
We have performed various fits of the curve depicted in fig. 4. The best one involves a 
combination of two terms of the form A/{r — t + c) n , with A, c, n determined by the 
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r-t 

Figure 4. The derivative of the function g(r — t) appearing in the asymptotic form of the field (eq. 
(5.1)). The model is the same as in fig. 1. 

data. Even though the data do not provide a unique answer, we find that a combination 
of powers with n ~ 3 — 5 gives the best fit. This behavior is consistent with the analytical 
results of ref. [2] for the theory described by the Lagrangian density (2.1). 

The DBI model with 62 = — 1 and A < 0, with an initial condition given by eq. (2.3), 
displays the second of the problems we described at the end of the previous section. A real 
solution over the whole range of r cannot be found after a certain time during the very 
early stages of the evolution. The problem appears near the extrema of (f> r , where 1 + \(f% 
vanishes. It is important to note that this term multiplies 4>u if the equation of motion is 
written in the form (3.3). In the following section we investigate the origin of this behavior. 

We next turn to the quartic model described by the Lagrangian (2.1). The most regular 
behavior is obtained for the model with <5i = — 1 and, therefore, A < 0. The evolution of 
the initial wavepacket is similar to that shown in fig. 1. In figs. 5, 6 we depict the solution 
at some early stage and at the latest time at which the solution is hyperbolic over the whole 
range of r (A > everywhere). The form of the solution is very similar to that depicted 
in figs. 2, 3 for the DBI model. In contrast to the DBI model with A > 0, the evolution at 
later times takes the system into the region where the equation of motion becomes elliptic 
(A < 0) within a certain r-interval. As we explained in the previous section, the boundary 
conditions provided by the scattering problem are not appropriate for the determination of 
a unique solution of the equation of motion when it becomes of mixed type (elliptic within 
a certain region, while staying hyperbolic in other regions). We find that instabilities tend 
to develop during the numerical integration when A becomes negative. The accuracy of the 
integration remains high during an initial stage at which the equation of motion is of mixed 
type. It diminishes gradually, until the numerical solution becomes unreliable at some time 
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Figure 5. The derivatives <fit (dashed) and <fi r (solid) of the nonlinear field, and the discriminant 
A of eq. (2.15) (solid grey), at two different times, in the context of the quartic theory, described 
by eq. (2.1), with 5i = —1, = 1. The initial wavepacket has A = 20, a = 1. The vertical dashed 
line denotes the classicalization radius. 

before the wavepacket reaches the classicalization radius (denoted by the vertical dashed 
line in figs. 5,6). 

Despite the above limitations, some useful information can be extracted from our 
solution. It is obvious from fig. 5 that the evolution of the wavepacket is similar to that in 
the DBI case with A > 0. The rear of the wavepacket becomes steeper and a shock front 
starts developing. This is indicated by 4> r becoming more negative. In the same time, a 
long tail is clearly visible at large values of r in fig. 6. The solid line in fig. 5 represents 
the discriminant A given by eq. (2.15). We have depicted the form of the wavepacket at 
the time when A first reaches zero at a certain value of r. 

Finally, the quartic model with A > displays the same behavior as the DBI model 
with A < 0. A real solution ceases to exist at some early stage in the evolution, at the time 
when the coefficient 1 — 3A<^ + Xtp^ of 4>tt i n the equation of motion (2.5) vanishes. 

6 Discussion and conclusions 

Despite some inherent problematic issues, our analysis of the idealized scattering process 
with exact spherical symmetry has provided concrete evidence for certain features of the 
classicalization scenario: 

• The classicalization radius of eq. (2.4) sets the scale for the onset of significant 
deformations of a collapsing classical configuration with large energy concentration 
in a central region. This behavior is consistent with the expectations of refs. [1-4]. 
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Figure 6. The nonlinear field cj) (solid) and the product 47rr 2 p, with p the energy density (dashed). 
The model is the same as in fig. 5. The vertical dashed line denotes the classicalization radius. The 
energy density is multiplied by 5 x 1CP 4 . 

• The equation of motion is a quasilinear partial differential equation of hyperbolic type 
at early times. At distances comparable to the classicalization radius, the nonlinear- 
ities become significant and can change the equation type. We discussed in detail 
an example in which the equation becomes essentially parabolic in certain regions of 
space. This is consistent with the analysis of ref. [6]. We also discussed an example 
in which the evolution makes the equation elliptic within a certain region. However, 
a real solution cannot be found at late times within the elliptic region. It seems 
likely that the initial conditions of the scattering process are not appropriate for the 
solution of the mixed type equation. 

• Shock fronts develop during the scattering process at distances comparable to the 
classicalization radius, consistently with the expectations of ref. [6]. 

• The most important observable feature of the classicalization process is the creation 
of an outgoing field configuration that extends far beyond the classicalization radius. 
This feature develops before the deformed wavepacket reaches distances of the order 
of the fundamental scale L*. 

Our study has also revealed several novel issues related to classicalization: 

• Our numerical analysis in the context of the DBI model indicates that the collapsing 
wavepacket can approach distances of the order of the fundamental scale before 
strong scattering appears. However, this could be a special feature of the DBI model. 
As we have seen, the equation of motion cannot become elliptic in this system, so 
that free propagation of wavelike configurations is not eliminated. 
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• The scattering resulting from the early stages of the classicalization process seems to 
be minimal within our analysis. The tail of the field configuration carries a negligible 
amount of energy because the corresponding modes are extremely soft. Within the 
DBI model, the bulk of the energy stored in the initial wavepacket can end up within 
a region of space with length scale comparable to L*. Similarly to the previous point, 
this behavior may be a special feature of the DBI model. 

• Our analysis does not provide evidence for the creation and subsequent decay of a 
quasistatic configuration that could be identified with the classicalon of refs. [1-4]. 
Instead, the classicalization scenario seems to be a fully dynamical process, whose 
defining feature is that it is characterized by a scale different from the naive length 
scale deduced from the higher-derivative terms in the Lagrangian density. A curious 
fact is that the static classicalons seems to exist in the two cases (quartic model with 
A > and DBI model with A < 0) in which a dynamical solution ceases to exist at 
an early stage of the evolution. 

• The most unexpected issue that emerged in our analysis is that the scattering problem 
may not have real solutions over the whole space within the context of the higher- 
derivative theories that support classicalization. 

In view of the above conclusions, a few more comments are in order concerning the 
accuracy of our analysis and the nature of the solutions. We have found that, in several of 
the versions of the equation of motion that we studied, an initial configuration representing 
a collapsing spherical wavepacket cannot evolve beyond a certain time. A real solution of 
the equation of motion ceases to exist for the initial condition of eqs. (2. 3), (2. 9). It is always 
possible that the inability to find a solution may be caused by the scheme we employed for 
the numerical solution of the equation of motion. For this reason, apart from the explicit 
leap-frog method we described in section 4, we reproduced our results through a second 
integration method, based on an implicit Crank-Nicolson scheme. The agreement is very 
good in the regions in which a solution exists, while both methods indicate the absence of 
a real solution in the same ranges of t and r. 

We would like to have a better understanding of the two cases in which a real solution 
ceases to exist at some early stage in the evolution, at the time when the coefficient of (fiu 
in the equation of motion (2.5) vanishes. These are the quartic model with A > and the 
DBI model with A < 0. It must be emphasized that this behavior is not associated with 
the presence of square-root singularities, similar to the ones appearing in the analytical 
solutions discussed at the ends of sections 2 and 3. The values of (fit, (fir are finite at the 
time at which a real solution ceases to exist. We can develop intuition on the problem by 
considering simple toy models that display similar behavior and can be solved explicitly. 
The simplest one is a harmonic oscillator of unit mass and non-standard kinetic terms. 

The analogue of the quartic model we discussed in section 2 is described by a La- 
grangian 
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with 5 = ±1 and £ > 0. The equation of motion is 



(1 -3<5£x 2 )x + w 2 x = 0. 



(6.2) 



Its first integral is the conserved energy 




(6.3) 



where we have assumed that initially the particle is at rest at x = xq. For this initial 
condition, the above equation gives 



For 5 = — 1 a solution exists at all times, describing a particle oscillating around the 
minimum of the potential at x = 0. For 5 = 1 a real oscillating solution exists at all times 
if Xq < (6£w 2 ) -1 . However, for Xq > (6£oj 2 )" 1 a real solution does not exist below the 
point with (6£ uj 2 ) . This is reached at the time when the coefficient of 'x in eq. 

(6.2) vanishes. We observed very similar behavior in the solutions described in section 5. 

The crucial property of the Lagrangian (6.1) with 5 = 1 is that the corresponding 
energy, given by eq. (6.3), is bounded from above when considered as a function of x. 1 
If the evolution starts at a suffciently large value of x, the initial potential energy can be 
larger than the maximal possible kinetic energy. As a result, the conservation of energy 
cannot be satisfied at some point during the evolution of the particle towards the minimum 
of the potential, and a real solution ceases to exist. The maximal value of the kinetic energy 
is obtained for dE/dx = x(l — 35^x 2 ) = 0. This explains why a real solution does not exist 
beyond the time at which the coefficient of 'x in eq. (6.2) vanishes. It is also interesting 
to consider the momentum conjugate to x, given by p = x(l — <5£x 2 ). This expression has 
a maximum for x 2 = l/(3<5£), given by \/ r S^p m = 2/(3\/3). There is no continuous real 
function x = x(p), defined over the whole real axis for p, that behaves like x ~ p for p — )■ 0. 

The field theoretical model of eq. (2.1) with 5\ = 1 and A > has very similar features 
with the toy model. The energy density of eq. (2.8) has a maximum as a function of (f>t at 
the point at which the coefficient of cfttt m ec i- (2-5) vanishes. This occurs at the time at 
which 1 — 3X(p 2 + X(j) 2 becomes zero for a certain value of r and the numerical integration 
stops converging. The origin of the problem is now clear: The energy cannot be conserved 
beyond this point and a real solution does not exist any more. The momentum density 
conjugate to <fi, given by ir = <^(1 — eft 2 + (ft 2 ), also has a maximum as a function of (fit when 
1 — 3\(fi 2 + Xcfi 2 = 0, similarly to the toy model. 

The analogue of the DBI model of section 3 is described by a Lagrangian 



Q5^uj 2 (x 2 -x 2 ) 
36l 



(6.4) 



L = - — Jl-5ix 2 




(6.5) 



lr The fact that the energy (6.3) is unbounded from below for x — > oo is not relevant for the classical 
evolution that we consider: The system starts with x = and the classical evolution breaks down at the 
finite value of x that corresponds to the "top of the hill" of the kinetic energy. 



with 5 = ±1 and £ > 0. The equation of motion is 



x + (1 - 5£x 2 f ,2 oj 2 x = 0, (6.6) 

and the conserved energy 

ZTi 1 , 1 2 2 1 , 1 2 2 /c>7\ 

= = H — w x = — H — uj x n , 6.7 

where we have assumed that initially the particle is at rest at x = xq. The above equation 
can be rewritten as 

l = h^ 2 (x 2 -x 2 ). (6.8) 



y/l-5ix 2 2 

For 5 = 1 a solution exists at all times, describing a relativistic particle oscillating around 
the minimum of the potential at x = 0. We observed very similar behavior in the solution 
depicted in figs. 1-3, which exists at all times up to the formation of the shock front. 

For 5 = — 1 a real oscillating solution always exists if X 2 , < (2£cj 2 ) -1 . However, 
for Xq > (2£w 2 ) -1 a real solution ceases to exist when the particle reaches the point 
with x 2 = Xq — (2£cj 2 ) -1 , at which x diverges. The origin of the problem can again 
be traced to the fact that the kinetic energy in eq. (6.7) is bounded from above when 
considered as a function of x. The maximum occurs for x — > oo. The conjugate momentum 
p = x/y/l + £x 2 also has a maximal value, equal to obtained for x — > oo. The 

corresponding field theory is the DBI model of eq. (3.1) with 82 = — 1 and A < 0. The 
partial derivative with respect to 4>t of the energy density of eq. (3.7) vanishes at the point 
where 1 + \<j) 2 , the coefficient of (fttt m eq. (3.3), becomes zero. The partial derivative 
with respect to (f>t of the momentum density tt = <fit / yl — \(f) 2 + A</> 2 also vanishes at the 
same point. These properties are responsible for the absence of a real solution of eq. (3.3) 
beyond this point, as the classical evolution cannot conserve the total energy of the system. 

In summary, our analysis has confirmed the existence of nontrivial physical behavior 
at distances much larger than the fundamental length scale in theories that support classi- 
calization. It has also revealed some peculiar features of certain underlying theories, which 
may lead to the absence of real solutions of the classical equations of motion. It remains 
to be seen whether the breakdown of the classical evolution is a pathology of such theo- 
ries, or whether it signals the transition to a regime in which the role of quantum physics 
is dominant. The one-dimensional toy models we considered in this section may provide 
crucial intuition on this issue. 
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